Atomic spin sensitive dissipation on magnetic surfaces 
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We identify the mechanism of energy dissipation relevant to spin-sensitive nanomechanics includ- 
ing the recently introduced magnetic exchange force microscopy, where oscillating magnetic tips 
approach surface atomic spins. The tip-surface exchange couples spin and atom coordinates, lead- 
ing to a spin-phonon problem with Caldeira-Leggett type dissipation. In the overdamped regime, 
that can lead to a hysteretic flip of the local spin with a large spin-dependent dissipation, even down 
to the very low experimental tip oscillation frequencies, describing recent observations for Fe tips on 
NiO. A phase transition to an underdamped regime with dramatic drop of magnetic tip dissipation 
should in principle be possible by tuning tip-surface distance. 



In a recent intriguing magnetic exchange force mi- 
croscopy experiment an exquisite magnetic atomic 
force sensitivity was demonstrated for an atomically 
sharp Fe magnetic tip over the (001) surface of antifer- 
romagnetic NiO 0] . Besides showing a different force for 
the two oppositely polarized surface Ni atoms — well ex- 
plained by the Fe-Ni exchange available from electronic 
structure calculations — the results also show a sur- 
prisingly different mechanical dissipation, with a gigantic 
excess of order 15-20 meV per cycle in the antiparallel Fe- 
Ni spin configuration, as compared to the parallel one. 

There is no existing theory of spin-dependent tip dis- 
sipation that one could use to understand not just this 
result but magnetically and site sensitive dissipation phe- 
nomena in general. Here we propose to use the magnetic 
exchange force microscopy study as a starting point. We 
search for a mechanism that i) can yield a magnetic dis- 
sipation of very large magnitude, similar to exchange en- 
ergies, per cycle; ii) is sensitive to the spin direction, 
and stronger for (nearly) antiparallel spin than a parallel 
one; iii) works down to the lowest frequencies. Partic- 
ularly puzzling is in fact the contrast between a large 
dissipation magnitude and the very low tip oscillation 
frequency (wtip ^160 kHz). At such a low frequency, 
one could expect a nearly adiabatic response, with very 
little mechanical energy transferred from the tip to some 
low-frequency excitations such as magnons, or perhaps 
phonons. Antiferromagnetic magnons, the first obvious 
choice, are immediately ruled out since, owing to strong 
dipolar anisotropy, the antiferromagnetic spin-wave dis- 
persion of NiO has a bulk gap A ~ 1.5 meV ~ 0.36 THz 
[H, and one at least as large at the surface 0-0- As a 
result, the oscillatory perturbation exerted on the sur- 
face spin is completely adiabatic — u>ti P <C A by more 
than 6 orders of magnitude — and direct dissipation in 
the spin-wave channel vanishes. Other strictly magnetic 
dissipation mechanisms involving mesoscopic scale phe- 
nomena, such as domain wall motion [8| , also appear un- 
applicable to the atomic scale tip-sample magnetic in- 
teraction. For example, a tip-induced magnetic domain 



with oscillating boundaries could be invoked to account 
for a low-frequency magnetic dissipation, but the forma- 
tion of such local domains is energetically unlikely, given 
the localized nature of the tip perturbation: Excluding 
a role of tip stray fields, simple model estimates suggest 
that the spin deformation near a perturbed surface spin 
should decay just a few atomic spacings away from the 
tip edge. We are left with acoustic phonons, certainly 
never gapped, both in bulk and at the surface. Here 
we know however that acoustic dissipation of a localized 
surface oscillation vanishes in linear response theory as a 
high power of frequency - the lattice can follow essen- 
tially adiabatically and harmonically a sufficiently slow 
and weak external perturbation. A large magnetic dissi- 
pation mechanism via phonons should therefore involve 
phenomena far from linear response. In this Letter we de- 
scribe the mechanism which we believe is at work here, 
and show that the nonlinear response is related to the at- 
tainment of a strong coupling overdamped spin-phonon 
state very well known in other contexts, giving rise to a 
single-spin hysteresis. That also suggests that by tuning 
down the perturbation intensity, a phase transition could 
be crossed from the overdamped to the underdamped 
regime, with a loss of hysteresis and a dramatic drop of 
dissipation. Hopefully, the present approach may serve 
as a prototype for nanoscale magnetic dissipation. 

Consider an oscillating Fe tip over a surface Ni spin 
Si. All neighboring spins remain essentially unper- 
turbed, "protected" as they are by the spin gap A. 
The potential felt by an f Ni atom at a distance z 
below the Fe tip differs from that felt by a i Ni 
and one can define a spin exchange potential V cx (z) = 
V^(z) — V^(z) (assuming the Fe tip to be i polar- 
ized). V ex can be estimated to yield an exchange force 
/ cx = -dV cx /dz of ~ 0.3 nN when the tip edge is 
closer than 3 A to the surface Ni This force pro- 
duces a small displacement u z (i) of the Ni atom from 
its equilibrium position and will result in a potential of 
the form — f c *Sfu z (i) (neglecting an unimportant spin- 
independent term). In terms of phonon creation (a^ ) 
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and annihilation (flks) operators (ks being wave vector 
and polarization of the phonon mode), we thus obtain 
a coupling of the Ni-spin to the Ni acoustic phonons 
of the form H< 



spinphonons 
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(a ks + a^ ks ), 
e 2 (ks); e z (ks) being 
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the eigenvector of the ks phonon mode. The equilibrium 
physics of the spin is dictated by the small frequency be- 
havior of the crucially important spectral density [lQ, El 



w ks )|e z (ks)| 2 /w k;j 



Wks)|Ak 



cx\2 



SAfJV ^ks 

From the standard Debye form for 
the low-energy acoustic phonons in three-dimensions, we 
find that the small-w limit of J{uS) is precisely Ohmic 
J(u) = h 2 auj + ■■■ with a = (/ cx ) 2 8M 3 fc ft 3 T » , where T D 
is the Debye temperature. An estimate, with / cx ~ 0.3 
nN, gives a value of a close to 1, which can be easily 
made > 1 by a slightly larger f CK or by a better ac- 
count of the (softer) surface phonon modes. The natural 
Ohmic behavior of this problem is a first important re- 
sult, since that i s by far the most interesting case, studied 
for decades 1(| 11 1, and it has been previously shown to 
arise in tip-surface interactions 12j . Note that the Ohmic 
coupling a depends on the square of the exchange force 
/ GX and is therefore dependent on the tip-atom distance 
z. We are thus led to the physics of a single spin — the 
surface Ni over which the tip is oscillating — in its proto- 
typical form, that of a driven Caldeira-Leggett (or spin- 
boson) Ohmic model. The model is known to possess two 
regimes, one underdamped and one overdamped, sepa- 
rated by a phase transition. In the underdamped regime 
the spin motion is relatively unaffected by the bath, and 
dissipation is small. In the overdamped regime — at- 
tained at a > 1 [lH — the spin is effectively "trapped" 
by the bath as schematically portrayed in Fig. [1] 

The tip-induced spin-flip processes involve in the over- 
damped case a new time scale r j~ 1 that can be much 
longer than the external driving period 2ttuj^. Once 
7 <C w-tip, magnetic dissipation will arise from a sort of 
single-spin hysteresis, similar to effects known in bistable 
models [T3 |. The overdamped model predicts three cru- 





Figure 1: Effective potential felt by the spin under the ef- 
fect of the bosonic bath for tip over spin-up and spin-down 
configuration. 



cial results regarding spin-dependent dissipation. First, 
dissipation is quantitatively large, because the order of 
magnitude of the hysteresis loop area is generally set by 
the tip-surface magnetic exchange, itself a large energy 
scale ~ 50 meV [3j. Second, dissipation will be strong 
when the magnetic tip is over a surface atom with an- 
tiparallel spin (left-hand side of Fig. [1]), and negligible 
over one with parallel spin (right-hand side of Fig. [TJ, 
because no tip-induced spin-flip is provoked in the lat- 
ter. Third, hysteretic dissipation should depend rela- 
tively weakly on 0Jf lp , and disappear only when the tip 
frequency is lowered below some very low frequency ~ 7, 
itself temperature-dependent. 

To describe the action of the tip on a given spin, we 
consider the driven spin-boson model (H = 1): 



H = % 



X 



(i) 



where a z and a x are the Pauli matrices, b\ and b v are cre- 
ation and annihilation operators for the mode v of the 
phononic bath of frequency w v , X — ^2 u ^u(b v + bl) 
is the bath operator to which a z is coupled, \ v being 
the previously introduced couplings, such that J(uj) = 



E„% - w„)|A„| 



(The slight time de- 



pendence of a during the tip oscillation is neglected.) 
In the NiO surface, £0 represents nearest-neighbor ex- 
change (positive and large) plus all sources of magnetic 
anisotropy (including dipolar contributions) , and z is the 
unperturbed direction of the surface atom spin due to 
dipolar anisotropy (a (211) direction). The external mag- 
netic tip acts in the xz plane at an angle 9 off the z 
axis, 

e(i) = ft- cos 6* sin 2 (o;tipt), A(t) — h sin 9 sin 2 (wti p t) . (2) 

We will use 9 and 9 + ir to describe Ni spins of oppo- 
site direction. The last term describes the free phonons. 
To avoid problems with divergences, we have as usual a 
high-frequency cutoff uj c for the spectral density, J(w) = 
awe" u/Ut , with u c > h, k B T. The details of the so- 
lution of this model, summarized below for the reader's 
convenience, are standard. The time evolution of the sys- 
tem is described by a standard real-time path-integral 
approach [l(| , taking care of the bath degrees of freedom 
through the influence-functional method and apply- 
ing the so-called noninteracting-blip approximation, valid 
in the a > 1 overdamped case for the observables of 
Actually, the same overdamped be- 
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havior could be achieved in the a < 1 regime, but only 
for certain values of perturbation and temperature. In 
terms of the free correlation function of the X opera- 
tor g(r) = S(r) + iR(r), with S(r) = f"du ^(1 - 

cos lot) coth ( 4p ) and R{t) = J^°duj ^^sinwr (where 
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P = we can compute the quantities Fo(t) = 

A 2 (t)j™dT e~ s ^ sin [R{t)\ sin [e(t)r] and G Q {t) = 
A 2 (t) / °°dr e~ s(T) cos [-R(r)] cos [e(*)r]. From these, as 
shown in [lti ] , in the overdamped regime a > 1 and with a 
low-frequency driving, the z component of the spin obeys 
a simple rate equation 



dt 



Mt)) = -Go(t) Mt))+F (t) , 



(3) 



that can be easily integrated. Applying the same proce- 
dure to (a x (t)), we get similarly 
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Fo(t)Go(t) 
A(i) 



(4) 



where Go(i) is defined as Go(t), but with sin[£(i)r] in 
place of cos [e(i) T ]- 

To uncover the new time scale, consider, e.g., the 
abrupt switching on of a perturbation at t = 0, with 
a constant value e(t) = e — £o and A(t) = A there- 
after. In this case Eq. (j3j) describes an exponential 
relaxation towards the equilibrium value (er z (oo)) = 
—F /Gq = — tanh(i/3e) [ll| with a decay rate 7 given, 
for low temperatures, in terms of the T function: 7 = 
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In the 
can take 



overdamped regime the relevant time scale 7 
large values, mainly due to the large cutoff frequency oj c 
being raised to a large exponent. Even under a compa- 
rably slow external perturbation the system can be out 
of its instantaneous equilibrium; and that is the origin of 
the hysteretic behavior. 

Consider now the two-level system of Eq. (JTJ under 
the external perturbation in ((2J. Figure [2] shows the 
time evolution of (<r z ) and (a x ) for a system in the ini- 
tial state {<7 Z ) = — 1, as obtained by numerical integra- 
tion of Eqs. © and (JU using an adaptive Runge-Kutta 
algorithm [17]. The system can clearly respond to the 
perturbation in its own time scale, yielding two different 
states when the perturbation is increasing and decreas- 
ing, a clear hysteretic behavior. The values of (a x (t)) are 
orders of magnitude smaller than those of (a z (t)) due to 
the large uj c . [This result is related to the universal be- 
havior of (a z (t)) with respect to lj c as opposed to the 
nonuniversality of (a x (t)) [la].] The dissipated energy 
per cycle is 



W = 



2-rr/u 



dt 
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dt 
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dA{t) 
dt 



, (5) 



which is the area of the hysteresis cycle in a (state- 
perturbation) (a z (t))-e(t) diagram. Figure [3] shows the 
hysteresis cycle for the z component of perturbations 
with different angles 6, together with the angular de- 
pendence of the hysteresis area (inset). When the tip 
and atom spins are (even roughly) opposite, the action 
of exchange to overturn the spin leads to a hysteresis loop 
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Figure 2: (Color online) Time evolution of (a z ) (black full line, 
left-hand axis) and (a x ) (dashed red line, right-hand axis) for 
£0 = 1, « c = 20eo, ot = 2.1, h = 1.5eo, # — 0.6, w t ip = 10 _8 £o, 
/3eo = 20. The dotted line shows the shape of the external 
perturbation. 



and a large dissipation; when they are nearly parallel, the 
loop collapses and correspondingly the magnetic tip dis- 
sipation drops. (Data for the x component, not shown, 
are negligible.) 

We may finally address the frequency and temperature 
dependence of the total magnetic dissipation. Loop areas 
(in steady state) for different temperatures as a function 
of frequency are shown in Fig. Hta). There clearly is an 
optimal frequency attaining maximal area. At excessive 
tip frequencies the spin remains effectively frozen in its 
trapped state; at very low frequencies the spin has plenty 
of time to relax and follow adiabatically the equilibrium 
value demanded by the tip: in both cases the loop area 
collapses. The inset of Fig. Ufa) shows how the opti- 
mal frequency increases with increasing temperature, re- 
flecting the T dependence of 7 shown in the forcing-free 
case. These results are consistent with what is known in 
the context of quantum stochastic resonance . Figure 
|3{b) contains the same data in the form of P = ojW, the 
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Figure 3: (Color online) Hysteresis cycle for the z component 
of the external perturbation for different angles 9 (see legend) 
for £0 = 1, w c = 100£o, a — 2.1, h — 1.3£o, Wtip = 10~ 12 £o- 
Inset: Angular dependence of the hysteresis area W. 
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Figure 4: (Color online) (a) Area of the hysteresis loop as 
a function of frequency for various values of /3 (see legend). 
Inset: Frequency of maximum area of the hysteresis loop as 
a function of the temperature T. (b) Dissipated energy as a 
function of frequency for various values of ft (see legend). 



dissipated power. At low frequencies P increases as a 
power law (roughly w 2 ), eventually reaching a plateau 
where the dissipation levels off over a wide frequency 
range. Other mechanisms will of course play a role at 
higher frequencies, but hysteretic dissipation is the only 
relevant nonlinear one that survives down to experimen- 
tally relevant low frequencies. 

The magnetic dissipation per cycle produced by the 
mechanism identified satisfies all desired requisites, since 
it is (i) large, and of the same order of magnitude of the 
antiferromagnetic exchange J k, 15-20 meV, (ii) vastly 
different for "up" and "down" Ni spins (assuming 9 = 
35A°, we get between the two a factor of 3 x 10 -4 ), and 
(iii) effective down to very low frequencies, ~ 7. Coming 
to the Fe-NiO data, we can now attribute the experi- 
mental dissipation of about 35 meV of the 4,-polarized 
tip oscillating over a 4. Ni spin to nonmagnetic mecha- 
nisms, that of about 50 meV over an t Ni in terms of 
the same mechanism plus a hysteretic magnetic dissipa- 
tion W w 15 meV, implying that the tip-surface coupling 
resulted in a > 1. 

The strong dependence of that coupling on / ox allows 
in principle for a reduction of a and a phase transi- 
tion from overdamped to undcrdampcd. In that case 
we would expect the faster evolution timescale to sup- 
press the hysteretic behavior, dramatically reducing the 



magnetic dissipation. 

In conclusion, our main novelties are that in sur- 
face magnetic tip dissipation problems, energy dissipa- 
tion should be mostly mechanical and non-spin-wave, 
since spin waves are generally gapped by anisotropy; that 
spin-dependent coupling to surface atomic motion and to 
phonons can lead to a sort of single-site magnetic hystere- 
sis; and finally that due to hysteresis the magnetic tip dis- 
sipation per cycle can be as large as intrasurface exchange 
coupling, as is seen experimentally. Further experimen- 
tal possibilities will be to test the frequency, temperature 
and 6 angle dependencies. Last but not least, the mod- 
ification of the coupling parameter with the tip-sample 
interaction force should in principle cause a phase tran- 
sition from an overdamped to an underdamped regime, 
with a strong suppression of dissipation at large distances 
and weaker couplings. We believe that these concepts 
should be of more general impact beyond the simple case 
treated here. 
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